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Generating  Functional  for  Electron-Phonon  Systems  from  Functional 
Integrals  with  One-Sided  Boundary  Conditions 

Joao  Florencio  Jr. 

Department  of  Physics  and  Astronomy,  University  of  Georgia,  Athens,  GA  30602 

We  derive  the  generating  functional  of  thermal  Green  functions 
for  an  interacting  electron-phonon  system  through  the  use  of  functional 
integrals  in  holomorphic  representations  (c-numbers  for  phonons  and 
Grassmann  numbers  for  electrons).  These  integrals  are  over  paths  obeying 
one-sided  boundary  conditions  at  either  0  or  0=l/kT,  and  do  not  require 
normalization  factors  unlike  previous  derivations.  We  recover  the 
Hatsubara  perturbation  formalism  without  using  any  results  from  the 
standard  operator  method. 


The  Matsubara  formalism  of  perturbation  theory  at  finite  temperatures 

plays  a  distinctive  role  in  the  theory  of  many-body  systems.1  There  one 

uses  operator  techniques  to  obtain  thermal  Green  functions  of  the  system 

as  a  perturbation  expansion.  Feynman  rules  can  then  be  used  to  evaluate 

each  term  of  the  expansion.  From  the  Green  functions  one  obtain 

correlation  functions,  which  are  closely  related  to  measurable  quantities. 

An  alternative  way  to  obtain  the  perturbation  expansion  is  by  the 

use  of  appropriate  generating  functionals  for  the  thermal  Green  functions. 

Their  forms  are  simple  enough  so  that  they  can  be  postulated  without 

further  justification  and  yet  they  are  known  to  reproduce  the  same  results 

obtained  by  the  operator  method.  There  are,  however,  some  attempts  to 

derive  generating  functionals  by  using  functional  integral  methods,  namely 

2  3 

in  the  work  of  Popov  and  Casher  et  al  .  These  works,  in  turn,  display 

2 

some  unsatisfactory  features.  In  the  work  of  Popov  the  field  operators 
are  required  form  the  beginning  to  be  periodic  (boson  case)  or 
antiperiodic  (fermion  case)  in  the  inverse  temperature.  There  requirements 
are  by  no  means  necessary  to  obtain  the  correct  Green  functions  --which 
are  indeed  periodic  (or  antiperiodic)  in  the  inverse  temperature--  as  one 

4 

can  easily  verify  this  by  using  the  standard  operator  method.  The 

3 

derivation  by  Casher  et  al  makes  use  of  an  incorrect  form  of  functional 
integrals  in  the  holomorphic  representation  for  bosons,  in  which  the 
integrals  used  are  over  periodic  paths.  Finally,  it  should  be  pointed  out 
that  none  of  these  works  give  the  explicit  expression  for  the  generating 
functional  of  the  ilatsubara  Green  functions  which  is  the  object  of  the 


present  work. 


In  this  paper  we  use  the  functional  integral  method  described  in 
Reference  6  to  obtain  the  generating  functional  for  the  thermal  Green 
functions  of  an  electron-phonon  system.  The  integration  variables  are 


c-numbers  for  the  phonons  and  Grassmann  numbers  for  the  electrons,  these 
variables  satisfying  one-sided  boundary  conditions  at  either  0  or  6. 

The  functional  integrals  used  here  are  known  to  reproduce  the  correct 
partition  functions  of  non-interacting  systems  of  bosons  and  fermions. 
They  were  also  used  to  obtain  non-perturbati ve  results  for  the  molecular 
polaron  with  quadratic  coupling  model. ^ 

Consider  a  system  of  electrons  and  phonons  in  a  homogeneous  medium 
of  volume  V.  For  simplicity  we  consider  electrons  with  the  same  spin 
and  phonons  belonging  to  a  single  branch.  The  electron  field  operators 
f(x)  and  ?+(x)  and  the  phonon  field  operator  $(x)  can  be  expressed  by 

4 

the  following  Fourier  expansions  : 


?(*)  ■  £  eikx  a.  , 

A  k  * 


(la) 


i\x)  =  I  e"’kx  a+  , 
/V  k  K 


and 


3(x)  -  -i  z fa.  (v^x  ♦  by1**) 

A  q  c  q  q 


(lb) 


(lc) 


The  operators  a^  and  a^  are  Fermi  annihilation  and  creation  operators, 

a  a  ^ 

respectively,  whereas  and  are  the  corresponding  Bose  operators. 
The  quantities  are  the  unperturbed  phonon  energies. 


3 


The  Hamiltonian  of  the  system  is  given  by 

H  =  H0  +  V  ,  (2) 

where 

"o  =  ^kVk  +  £“q(Sqeq  +  H)  (3) 

K  q 

is  the  non-interacting  part,  and 

V  =  Jdx  v(v+(x),v(x),£(x))  (4) 

is  the  interaction  between  electrons  and  phonons.  In  the  above,  ek  are 

A 

the  unperturbed  single  electron  energies,  and  v  is  the  interaction 
Hamiltonian  density. 

Now  we  introduce  the  temperature-dependent  external  sources 
£,  £*  and  n,  which  interact  linearly  with  the  electron  and  phonon 
operators  as  follows: 

H j ( t )  =  -  Jdx  (v+(x)5(x,t)  +  £  (x,th(x)  +  o(x)n(x,x))  .  (5) 

The  sources  £  and  £*  anti  commute  with  the  electron  operators  and  among 
themselves,  whereas  n  commutes  with  the  phonon  operators.  These  sources 
need  not  have  actual  physical  meaning,  and  at  the  end  we  let  them  to 
vanish. 

The  grand-Hamiltonian  K  of  the  system  in  the  presence  of  the  sources 
is  given  by 

K( t)  =  H  +  Hj(t)  -  uN  , 


(6) 


4 


where  u  is  the  electron  chemical  potential,  and  N 


k  is  the 


total  number  of  electrons  operator. 

The  generating  functional  is  defined  as  the  following  functional 
of  the  external  sources: 


Z[5*,5.n]  =  Tr  p[c*.C,n]  , 


where 


p[5*.5tn] 


6. 

T  exp{-  /  K(x)dT} 
T  0 


(7) 


(8) 


is  the  density  matrix  in  the  presence  of  the  sources.  Here,  is  the 
temperature  ordering  operator.  The  trace  operation  can  be  performed  in 
the  holomorphic  representation  for  the  electrons  and  phonons,  in  which 
the  electrons  are  represented  by  Grassmann  numbers  --a^  and  a£--  and 
the  phonons  by  c-numbers  --bq  and  b*.  For  details  of  holomorphic 

O 

representations  see,  e.g.,  Faddeev  and  Slavnov  and  references  therein. 

In  the  holomorphic  representation,  the  right  hand  side  of  Eq.(7) 
is  given  by 

db*db 

Tr  o[s*,5,n]  =  //  jn(b*,a*,b,a)  nexp(-b*bq)  ■  3  q 

x  n  exp(-aka*)  dakda*  ,  (9) 

k 

where  o,*  ,  (b*,a*,b,a)  is  the  kernel  of  the  density  matrix  operator, 

and  the  subscripts  denote  its  functional  dependence  on  the  external  sources. 

The  density  matrix  kernel  is  represented  by  the  following  functional 

,6 

integral  : 


p r *  r  _(b*,a*,b,a) 

s  t n 


/b$<B)-b$  exp{Eb*(B)bq(6>  *  Ia;;(e)<>k(e) 

bq(°)=b, 


6 

+  /  (  -  zb*(u)b  (u)  -  za*(u)a.(u)  -  K(b*(u) ,a*(u) ,b(u) ,a(u) ,u))du} 
0  q  q  q  k  K  K 


db*(u)db„(u) 

x  (  n  —9 - 9 — .)  (  n  da£(u)dak  (u))  ,  (10) 

q,u  2iri  k,u 

8  '** 

where  K  is  the  normal  symbol  of  the  operator  K,  that  is, 

K  »  t(ek-M)a*ak  +  Zwq(b*bq+Js)  +  Jdx  v(t*,¥.0  -  Jdx(t*<;  +  £*Y  +  *n). 
k  q 

(ID 


In  the  above  expression,  v  is  the  normal  symbol  of  the  interaction 
density  v,  and  y*,  y  and  <t>  are  the  normal  symbols  of  the  field  operators 

A  A 

y  ,  Y  and  ♦,  respectively.  Notice  that  the  functional  integral  is  neither 
periodic  nor  anti  periodic  in  the  inverse  temperature,  rather  it  obeys 
one-sided  boundary  conditions  at  either  0  or  6. 

The  interaction  part  can  be  removed  from  the  integrand,  by  using 
the  following  identity: 
s 

exp{  /  du/dx(-v(Y*,Y,$)  +  Y*(x,u U(x,u)  +  s*(x,u) y(x,u)  +  4>(x,u)n(x,u))  } 
0 

s  s 

=  exp{  -  /  du/dx  >  exp{^/  du/dx(Y*(x,uk(x,u)+£*(x,u)Y(x,u) 

+$(x,u)n(x,u)) }  ,  (12) 


where  we  have  made  the  following  replacements  in  v. 


(right  derivative). 


and 


¥* 


0 

6C 


(left  derivative) , 


The  exponential  form  containing  v  on  the  right  hand  side  of  Eq.(12) 
is  to  be  understood  as  a  power  series  expansion. 

By  using  Eq.(12),  the  functional  integrals  for  the  electrons  and 
phonons  become  decoupled  from  each  other.  The  density  matrix  kernel 
assumes  the  following  form: 


6 


0t*.t,n(b*'a*>b'a)  =  du'dx  V(S?'S?*^)}  °n 


ph 


(b*»b) 


(13) 


where  pph  and  are,  respectively,  the  density  matrix  of 

n  s  » s 

non-interacting  phonons  and  electrons  in  the  presence  of  external  sources. 
They  are  given,  respectively,  by 


oPh(b*,b) 


=  W)=b*  6XP  E{b5(6)bq(6)  +  J  (-b^u)bq(u)  -  UJq(bq(u)bq(u)+^) 


q  '  q 

bq(0)=bq 


db*(u)db  (u) 

-  (b*(u)y  (u)  +  b  (u)y*(u)))du}  n  —3 - 9 - 


(14) 


7 


and 

el  6 

p®;^(a*,a)  =  /a*(6)=a*  exp  r { a*( e )afc ( e )  +  /  (-a*(u)ak(u)  -  (ek-p)a£(u)ak(u) 

k  k  k  0 

ak<°)=ak 

-  (a^(u) Ck(u)  +  s£(u)ak(u)))du}  n  da*(u)dak(u) ,  (15) 

K » u 


where 

Yq(u)  =  -  /dx  e“iqxn(x,u)  , 

y*(u)  =  -  /dx  eiqxn(x,u)  , 

ck(u)  = - -  /dx  e_lkxc(x,u) , 

*  AT 


Cf(u)  = - -  /dx  eikxc*(x,u)  . 

K  /V 


(16a) 

(16b) 

(16c) 


(16d) 


Both  of  the  integrals  (14)  and  (15)  are  Gaussian  and  can  be  evaluated 
exactly.  Their  solutions  are  simply  the  values  of  their  integrands 
at  the  paths  which  extremize  their  actions,  i.e.,  the  arguments  in 
the  exponential  forms®.  The  extremal  conditions  yield  the  equations 


b  (u)  +  wqbq(u)  +  Yq(u)  =  0, 

V0)=bq 

(17a) 

b*(u)  -  wqb*(u)  -  Y*(u)  =  0, 

(17b) 

8 


$ 

I 


m 


aR(u)  +  Uk-u)ak(u)  +  ;k(u)  =  0,  ak(0)=ak 


(17c) 


a£(u)  -  (ek-y)a*(u)  -  ?*(u)  =  0, 


Their  respective  solutions  are 


ak(B>=ak- 


(17d> 


b  (u)  =  b  e'V  -  /  e-“q(u's)Y  (s)ds, 

q  q  o  q 


b*( u )  =  b*e  q(u_6)  -  / Vq(u's)Y*(s)ds, 

a.  (u)  =  a,  e”^Ek"M^u  -  /  e~^k~u^u~s\.  (s)ds. 


a*(u)  =  a£e^k~w^u"6)  -  /  e^Ek‘u^u"s^^(s)ds. 


(18a) 


(18b) 


(18c) 


(18d) 


Notice  that  neither  these  extremal  paths  nor  their  respective 
equations  are  complex  conjugate  pairs. 

The  values  of  the  integrals  (14)  and  (15)  are  then  given,  respectively. 


.  S(J  -  „  3 

>^n(b*,b)  =  exp{-  E— }expE{b*bqe  6ujq  -  b*e  3"q  f  exqJYq(u)du 

-  bn  /  e~UJqUY*(u)du  +  /  du  /  ds  e'uq^u_s^Yq(u)Yq(s) }  (19) 


0  0 


Pc*>5(a*,a)  =  expE{a*ake_tne 
-  (  /  e“^£k"M^uc*(u)du)a, 

0  K  K 


k'v>  -  aJe-^V^  /  euk^V(u)du 
K  0  K 

+  /du  /  ds  e"(ek‘p)(u's)^(u);.(s)} 

0  0  K  K 


(20) 


These  are  the  final  results  for  the  density  matrix  kernels  of 
non-interacting  phonons  and  electrons  in  the  presence  of  external 
sources.  Next  we  use  the  trace  formulae  to  obtain  the  following 
partition  functions: 
n  nh  db*db 

ZphH  =  /<  (»*.b)exp(-Eb*b  >  n— 3— 3- 
K  q  M  H  q  2tt1 


q  3  u  /  \ 

=  (  n  *1.  )exp(s((N°+l)  /  du  /  ds  e‘Vu“s)Y*(u)Y 0(s) 
q  1-e  %  q  q  0  0  q  q 

8  8 

♦  /  du  /  ds  e'“q(u's),*(u)v(s)))  (21) 

q  0  u  q  q 


and 

ze-j  [€*•€]  =  /p®I>5(a*,a)exp{-zaka*}  ndakda* 

k  k 

=  (  n(l+e'3^E:k'u^))exp{zf  /  du  /  ds  e“^k"u^u“s^*(u);.  (s) 
k  k  0  0 

-  NF  /du  J3ds  e"(ck~u ) ^u_s > ;*(u) ck(s ))  }  , 

K  n  n  n  n 


(22) 


10 


n.f=  1 


k  e6lVu)  +  1 

are,  respectively,  the  familiar  Bose  and  Fermi  distribution  functions. 

Notice  that  the  dependence  of  the  partition  functions  on  o,  £*  and  £ 

is  implicitly  understood  in  the  above  expressions  due  to  Eqs.(16a-d). 

We  now  use  the  following  identities 

B  u  0  8 

KnJJ+U  J  du  /  ds  e’“q(u-s)r5(ii)rq(s)  »  6uJ  ds  e-“qlu-sl-,*(u)Yq(s), 


0  0 


ENq  /“du  JS ds  e'“q(u's)-,;(u)v  (s)  »  /“du  1° ds  e“q I u's I T*(u )y  (s) 


0  0 


(23a) 


8  8 


q  H  0  u 


q  '  q 


q  H  0  0 


’q'  ' 'q' 
(23b) 


and 


l  J  du  /  ds  e"^ek"y^u_s^A(u)c.  (s)  =  l  /  du  /  ds  e"^Ek"u^u~s^^(u) c.  (s)e(u-s 
kOO  **k00  *  * 

(23c) 

where  e  is  the  step  function,  together  with  Eqs.(16a-d)  to  express  the 
partition  functions  (21)  and  (22)  in  the  following  forms: 


3  6 


Zph[n]  =  (  n  )exP(|  /  du  /  ds/dx/dy  n(x,u)|  l^.f(N®+l)eiq(x'y)e‘uq  lu'S  I 

v  q  1-e  q  "0  0  q 


*  ^e'<l<*-y>e“,l“-5l)„(y.S)>  (24) 

and 

6  8 

Zp, [£*,£]  =  (n(l+e'6^£k'w^))exp{  /  du  /  ds/dx/dy 
ei  k  0  0 

x  £*(x,u)i  Ee1'k(x-y)e-(ck_lj)(u's)(e(u-s)-N!')£(y,s)}  .  (25) 

v  k  K 


Notice  that  the  unperturbed  thermal  Greem  functions  for  electrons  and 
phonons  appear  explicitly  in  the  above  expressions.  In  addition  , 
the  pre-exponential  forms  are  simply  the  unperturbed  partition  functions. 
More  explicitly, 

g  g 

Zph [n]  s  Z°h  exp{-  j  /  d uj  ds/dx/dy  n(x,u)D^0^ (x,u;y,s)n(y,s) }  (26) 
and 


Zel [«*•«! 


g  g 

zJi  exp{-  /  du  /  dsjdx/dy  c*(x,u)G^0^ (x,u;y,sK(y,s) }  , 
e  0  0 

(27) 


where 


2°  -  „ 

»h  q  l-e-s"q  ’ 


Z°,  -  "(1  *  )  , 


(28) 

(29) 


D^(x,u;y,s) 


s - ((N®+1)eiq(x”^)e~u,q ,u_s  1  +  N®ei<»(x‘y)euqlu‘sl)  , 

^  *  n  ■  '  • 


(30) 


and 


G^(x,u;y,s)  = 


1  reik(x-y)e-(ek-u)(u-s)(g(u_s)  _  f,F} 
v  k  K 


(31) 


The  above  expressions  for  the  unperturbed  phonon  Green  function  0^ 

and  the  electron  Green  function  G^  are  the  same  found  in  textbooks.^ 
Notice  that  they  were  obtained  here  strictly  from  functional  integral 

methods.  There  was  no  assumption  made  concerning  periodicity  or 

antiperiodicity  of  the  field  operators  in  the  imaginary  time,  yet  these 


Green  functions  do  obey  periodicity  — D^  —  and  anti  periodicity  — G^-- 
conditions. 

Finally,  by  combining  the  above  results  we  arrive  at  the  final 
expression  for  the  generating  functional  for  the  interacting  system: 

Z[5‘,S,n]  ’  Z°  exp{-  Q/SdU/dx 

1  S  6  /  Q  \ 

x  exp{-  7  /  du  /  ds/dx/dy  n(x,u)Dv  ‘ (x,u ;y ,s)n(y ,s) 

^00 

0  0 

-  /  du  /  ds/dx/dy  $*(x,u)G^(x,u;y,s)e(y,s) }  ,  (32) 

0  0 

where  Z^ 

wnere  L  Sh  eT 

The  partition  function  is  obtained  by  first  expanding  the 
form  containing  the  interaction  v  in  a  power  series,  and  then 
performing  the  functional  derivatives  required  in  each  term  of  the  series 
and,  finally,  setting  the  sources  equal  to  zero.  Of  course,  one  still 
is  left  with  integrals  over  the  spatial  coordinates  and  the  inverse 
temperature.  This  procedure  amounts  to  summing  all  possible  diagrams 
with  no  external  lines. 

The  two-point  Green  functions  for  the  electrons  and  phonons 
are  obtained,  respectively,  by 


G( Xj  ,  T  ^  1X2 , 


t2 


)  =  - 


1  6zZ 


5*=5=n=0 


(33) 


and 


0( X ^ , t ^ j X2  ,  t £ )  ” 


I  6zZ 

Z  6n(x^,Tj)5n(x2»T2) 


C*=C=n=0 


(34) 


V*. 
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The  diagram  expansions  of  these  Green  functions  will  contain  only 
connected  diagrams  with  two  external  lines.  As  usual,  each  line 
in  a  given  diagram  is  associated  with  an  unperturbed  Green  function 
and  the  vertices  are  associated  with  the  interaction  v.  The  coordinates 
x1  and  X£  denote  the  free  ends  of  the  external  lines  at  and  t2» 
respectively.  In  these  expansions,  Wick's  theorem  is  automatically 
satisfied.  Finally,  it  should  be  pointed  out  that  Green  functions 
with  many-points  can  also  be  obtained  from  (32). 

To  summarize,  we  have  derived  the  explicit  expression  of  the 
generating  functional  for  the  temperature  dependent  perturbation 
theory  of  an  electron-phonon  system.  The  form  (32)  which  insofar 
has  only  been  used  heuri stical ly  is  now  indeed  justified.  One  should 
also  notice  that  no  use  of  prior  results  from  the  operator  method 
was  required  to  the  derivation  of  the  generating  functional.  In  a 
forthcoming  article  we  use  techniques  of  functional  integrals  in  the 
holomorphic  representation  to  show  how  to  obtain  non-perturbative 
expansions  for  the  partition  functions  of  systems  of  bosons  and  fermions. 

ACKNOWLEDGMENTS 

The  author  is  grateful  to  Professor  B.  Goodman  for  several  valuable 
discussions.  He  is  also  indebted  to  Professor  M.  H.  Lee  and  the  Office 
of  Naval  Research  for  the  financial  support. 


14 


REFERENCES 

1.  T.  Matsubara,  Prof.  Theor.  Phys.  14_,  351  (1953). 

2.  A.  Casher,  D.  Lurie,  and  M.  Revzen,  J.  Math.  Phys.  9_,  1312  (1968). 

3.  V.  N.  Popov,  Functional  Integrals  in  Quantum  Field  Theory  and 
Statistical  Physics  (Dordrecht:  Reidel,  1983). 

4.  See,  e.g. ,  A.  A.  Abrikosov,  L.  P.  Gorkov,  and  I.  E.  Dzyaloshinsky, 

Methods  of  Quantum  Field  Theory  in  Statistical  Physics  (Englewood  Cliffs: 
Prentice-Hall,  1963). 

5.  In  the  passage  from  Eq.(7)  to  Eq.(8)  of  Ref.  2,  the  authors  mislabel 
some  of  the  quantities  involved:  There  is  no  justification  for  setting 
<jg  =  in  some  terms  while  leaving  the  combination  ogc^  unchanged. 

6.  J.  Florencio  and  B.  Goodman,  Phys.  Rev.  B  3£,  3634  (1986). 

7.  J.  Florencio  and  B.  Goodman,  Phys.  Rev.  B  34,  3639  (1986). 

8.  L.  D.  Faddeev  and  A.  A.  Slavnov,  Gauge  Fields  (Reading:  Benjamin/ 
Cummings,  1980). 

9.  The  form  (30)  can  be  found,  for  instance,  in  Ref.  1,  p.  103,  whereas 
the  form  (31)  is  in  G.  Mahan,  Many-Particle  Physics  (Plenum:  New  York, 
1981),  p.  130. 


The  time  evolution  in  the  XY  model  is  expected  to  be  quite  different 

from  that  in  the  TI  model  since  the  two  models  are  inequivalent.  If,  however, 

the  strength  of  the  external  field  B  in  the  one-dimensional  TI  model  is 

made  to  approach  the  coupling  constant  J,  the  time  evolution  becomes 

identical  to  that  in  the  one-dimensional  XY  model.  We  have  shown  that 

when  B  s  J,  the  geometric  structure  of  the  Hilbert  spaces  of  S *  in  the 

0 

TI  and  XY  models  become  identical.  The  two  structures  are  not  identical 
at  any  other  values  of  B.  It  has  been  established  that  the  form  of  the 
time  evolution  is  a  unique  function  of  the  geometric  structure  of  the 
realized  Hilbert  spaced  Hence  one  can  understand  the  origin  of  the 
coincidence  in  the  time  evolution  of  Sx.  at  this  particular  value  of  B. 

J 

In  addition,  the  realized  recurrence  relation  turns  out  to  be  simple 
enough  to  permit  an  analytic  solution  for  the  entire  family  of  the 
relaxation  functions. 

For  the  XY  and  TI  models  in  higher  dimensions,  e.g.  d  =  2,  one  does 
not  presently  know  anything  about  their  time  evolution.  While  the 
analytic  solution  for  it  poses  a  challenge  to  any  method,  there  is  a  simpler, 
perhaps  more  immediate  problem.  That  is,  whether  the  coincidence  in  the 
time  evolution  is  unique  to  one  dimension,  or  whether  it  also  exists  in 
higher  dimensions  possibly  with  some  d-dependent  values  for  B.  This  problem 
can  be  readily  tackled  by  the  method  of  recurrence  relations,  since  the 
recurrants  can  be  calculated  in  higher  dinensions.  According  to  the  method 
of  recurrence  relations,  it  is  then  sufficient  to  compare  the  geometric 
structures  of  the  realized  Hilbert  spaces  and  to  find,  it  exists,  the 
value  of  B  for  which  the  two  structures  become  identical  to  each  other. 


